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Abstract: Wc consider the effective area Acfr of photonic crystal fibers 
(PCFs) with a triangular air-hole lattice in the cladding. It is first 
of all an important quantity in the context of non-linearities, but it 
also has connections to leakage loss, macro-bending loss, and numerical 
aperture. Single-mode versus multi-mode operation in PCFs can also 
be studied by comparing effective areas of the different modes. We 
report extensive numerical studies of PCFs with varying air hole size. 
Our results can be scaled to a given pitch and thus provide a general 
map of the effective area. We also use the concept of effective area to 
calculate the "phase" boundary between the regimes with single-mode 
and multi-mode operation. 
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1 Introduction 

Photonic crystal fibers (PCF) constitute a new class of optical fibers which has revealed 
many surprising phenomena and also has a big potential from an application point of 
view (for recent special issues on photonic crystal fibers, see e.g. Refs. [1, 2]). Here we 
consider the type of PCF first studied in Ref. [3, 4] which consists of pure silica with a 
cladding with air-holes of diameter d arranged in a triangular lattice with pitch A, see 
Fig. 1. For a review of the basic operation we refer to Ref. [5]. 

The effective area is a quantity of great importance. It was originally introduced as 
a measure of non-linearities; a low effective area gives a high density of power needed 
for non- linear effects to be significant [6]. However, the effective area can of course also 
be related to the spot-size w through A^ff = ttw'^, and thus it is also important in 
the context of confinement loss [10], micro-bending loss [11], macro-bending loss [12], 
splicing loss [13], and numerical aperture [14]. 

Strictly endlessly single-mode operation of PCFs [15] is possible for d < d* ~ 0.45A 
[5] , but even for larger air holes single- mode operation is possible for wavelengths above 
a certain cut-off A*. We demonstrate that the effective area of the second-order mode 
is a useful concept in determining this cut-off for a given hole size d. For d < d* we 
recover the endlessly single-mode regime with A* — *■ 0. 

The paper is organized as follows: In Sec. II we introduce the concept of the effective 
area and In Sec. HI we report numerical results for PCFs of the type shown in Fig. 1. 
Finally, discussion and conclusions are given. 



(1979). 




Fig. 1. Schematic of end-facet of an all-silica photonic crystal fiber. The microstruc- 
tured cladding consists of air holes of diameter d arranged in a triangular lattice 
with pitch A. The silica core guiding the light is formed by the "missing" air hole 
(indicated by the dashed circle). 



2 Eigen-modes and effective areas 

We consider PCFs realized in dielectric materials with no free charges or currents. From 
Maxwell's equations it then follows that the iJ-field is governed by the general wave 
equation (see e.g. Ref. [16]) 



V X 



e(r)V X H^{r) = {uj/cfH^{r) 



(1) 



where e is the dielectric function and lo is the frequency of the harmonic mode, Hi^{r, t) = 

The PCF has translation symmetry along the fiber axis (z-axis) and for e{r) = 
the solution is of the form 



±i/3„(w)z 



(2) 



where hn is the transverse part of the nth cigcn-modc and /3„ is the corresponding 
propagation constant at frequency u. Often one will specify the free-space wavelength 
A rather than the frequency uj = c(27r/A). 

The effective area associated with the nth eigen-mode is given by [6] 



^cff,„(A) = [ jdr^Ur^)\y j dr^I^ir^), (3) 

where /n('r±) = \hn{r±)\'^ is the intensity distribution. It is easy to show that for a 
Gaussian mode h{r±) oc e^'^'"^/™-' of width w the effective area is A^s = nw^. Applying 
Eq. (3) to closc-to-Gaussian modes in some way corresponds to a Gaussian fit averaged 
over all angular directions. In Fig. 2 we illustrate this by an example. For d/A = 0.3 
and A/A = 0.48 we compare the intensity / of the fundamental mode (n = 1 or 2) to 
the intensity Iq of the corresponding Gaussian with the width calculated from Aaff. As 
seen the over-all intensity distribution is reasonably described by the Gaussian with a 
width calculated from the effective area. 

The low-intensity deviations from a Gaussian intensity distribution are coursed by 
the six-fold symmetry of the air hole lattice; the field extends slightly into the six silica 
bridges formed by the six holes nearest to the core. Similar modal properties have been 
found experimentally, see e.g. Refs. [3, 7, 8]. 

Below we mention a few examples of some of the different phenomena and quantities 
which can be at least partly quantified from knowledge of the effective area; the non- 
linearity coefficient, the numerical aperture, the macro-bending loss coefficient, and 
splicing loss. In Sec. Ill we also discuss how a second-order mode cut-off can be identified 
using the concept of the effective area. 

2.1 Non-linearity coefficient 

The non-linearity coefficient 7 is given by [6] 



cAeS XAeff 

where n2 is the nonlinear-index coefficient in the nonlinear part of the refractive index, 
Sn = n2\E\'^ . Knowledge of A^s is thus an important starting point in the understanding 
of non-linear phenomena in PCFs. Due to the high index contrast between silica and 
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Fig. 2. Comparison of the real intensity I{r±) to the Gaussian intensity /G(rx) 
with a width obtained from A^ff. The upper left panel shows Ic/Iy the upper right 
panel shows /, the lower right panel shows Iq, and the lower left panel shows the 
dielectric function e with d/A = 0.3. The fields are calculated at A/A = 0.48. 

air the PCF technology offers the possibihty of a much tighter mode confinement (over 
a wide wavelength range) and thereby a lower effective area compared to standard-fiber 
technology. Furthermore, the microstructured cladding of the PCFs also allows for zero- 
dispersion engineering. For a recent demonstration of a highly nonlinear PCF with a 
zero-dispersion at A = 1.55 see Ref. [9]. 

2.2 Numerical aperture 

Also the numerical aperture (NA) relates to the effective area. For a Gaussian field 
of width w one has the standard approximate expression tan^^ ~ X/ttw for the half- 



divergence angle 6 of the light radiated from the end-facet of the fiber [17]. The corre- 
sponding numerical aperture can then be expressed as 

NA = sin6' ~ [l + TiA^al>?y^''^ ■ (5) 
In Ref. [14] this was used in a study of the numerical aperture of PCFs. 

2.3 Macro-bending loss 

For estimation of the macro-bending loss coefficient a the Sakai-Kimura formula [18] 
can be applied to PCFs. This involves an evaluation of the ratio A\/P, where Ae 
is the amplitude coefficient of the field in the cladding and P the power carried by 
the fundamental mode. The Gaussian approximation gives A^/P = \/Acs [19] and in 
Ref. [12] this was used to calculate the macro-bending loss in PCFs based on fuUy- 
vectorial eigenmodes of Maxwell's equations. 

2.4 Splicing loss 

Also splicing loss can be quantified in terms of the concepts of effective areas. The 
splicing of two aligned fibers with effective areas Aeff,i and ^off,2 will have a power 
transmission coefficient T < 1 given approximately by 

due to the mismatch of effective areas. 

3 Numerical results 

We solve Eq. (1) numerically and from the eigenmodes we calculate the corresponding 
effective area from Eq. (3). The fully- vectorial eigenmodes of Maxwell's equations are 
computed with periodic boundary conditions in a planewave basis [20] . This approach 
provides the eigenmodes on a discrete lattice and the integrals in Eq. (3) are then 
obtained by sums over lattice sites inside the super-cell. For the dielectric function we 
have used e = 1 for the air holes and e = (1.444)^ = 2.085 for the silica. Ignoring the 
frequency dependence of the latter the wave equation, Eq. (1) becomes scale-invariant 
[16] and all the results to be presented can thus be scaled to the desired value of A. 

3.1 Fundamental mode 

In Fig. 3 we show the effective area of the fundamental mode as a function of wavelength 
for different hole sizes. As expected the mode becomes more confined - lower effective 
area - for increasing air hole size. In general we find that the effective area of the 
fundamental mode is of the order of A^ with a prefactor which depends slightly on the 
air hole size. In fact, scaling the results by a factor d/A we find numerical evidence that 

Aoff oc (A/d) X A^ + 0{\/d), (7) 

and as demonstrated in the lower panel of Fig. 3 the prefactor is of order 0.5. 

Fig. 3 is a general "map" of the effective area which can be used in calculating 
the different quantities discussed in Sec. II as well as for designing single-mode PCFs 
with specified modal properties. If for example the value of d/A is given (which is 
often the case during stack-and-puU production of PCFs) a given effective area can be 
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Fig. 3. The upper panel shows the effective area as a function of wavelength for 
different hole sizes. The lower panel shows the same data as in the upper panel, but 
with the vertical axis scaled by the factor d/A and the horizontal axis scaled by the 
factor (d/A)-i. 

obtained by scaling the fiber structure, i.e. the pitch A. The Corning SMF28 standard 
fiber has an effective area ^off ~ 86 /J,m^ at A = 1550 nm and Fig. 3 suggests that a 
comparable effective area can be realized by the PCF technology using e.g. d/ A ^ 0.25 
and A 6.5 fim. 

3.2 Second-order mode 

In Fig. 4 we show a calculation of the effective area of the second-order mode (n = 3) for 
a PCF with d/ A ~ 0.5 air holes. For short wavelengths the effective area is of the order 
A^ with the mode confined to the core of the PCF. At long wavelengths the effective 
area diverges corresponding to a dclocalizcd cladding mode. Even though the transition 



to a cladding mode is not abrupt the high slope of the effective area allows for a rather 
accurate introduction of a cut-off wavelength A*. This cut-off wavelength is indicated 
by the crossing of the dashed line with the horizontal axis. 

3.3 Single-mode versus multi-mode operation 

Conditions for single-mode operation of PCFs turn out to be very different from stan- 
dard fiber technology and for not too big air holes, d < d* , they can even be endlessly 
single- mode [3] . The reason is that the effective cladding index has a wavelength depen- 
dence which prevents the core from guidance of higher-order modes while still guiding 
a fundamental mode. In Ref. [5] d* ~ 0.45A was suggested from standard-fiber consid- 
erations. 

For d > d* a second-order mode is guided for wavelengths shorter than a certain 
cut-off and above this cut-off the PCF is single-mode. This cut-off of the second-order 
mode can be found from studies of the corresponding effective area; below the cut-off 
the effective area Acs.s is finite and comparable to and above the cut-off the effective 
area diverges (in a super-cell calculation it approaches the area of the super-cell). 

Comparing Fig. 4 to Fig. 3 for d/A = 0.5 it is seen how we for wavelengths shorter 
than A* ~ 0.35A have a multi-mode PCF whereas the PCF is in the single-mode regime 
for A > A* . Carrying out the same analysis for different hole-sizes allows for constructing 
a "phase" diagram, see Fig. 5. The data points indicate calculated cut-off wavelengths 
and the solid line is a fit to the function 

f{x) = a{x - xof, xo = 0.45, (8) 

with the fitting coefficients a ~ 1.34 and P ~ 0.45. 

The solid line defines a "phase" boundary; above the line the PCF is single-mode 
and below it is multi-mode. Furthermore there is an endlessly single-mode regime for 
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Fig. 4. The effective area of the second-order mode for d/A = 0.5. For short wave- 
lengths Acji_3 ~ and the mode is guided in the core. For high wavelengths 
the mode becomes a cladding mode with effective area approaching the area of the 
super-cell used in the calculation. The crossing of the dashed line with the horizontal 
axis indicates the cut-off for the second order mode. 



d < d* ^ 0.45A [5]. In principle the effective area approach aUows for an independent 
determination of d* but it should be emphasized that the numerical efforts needed to 
resolve the cut-off increase dramatically when d approaches d* . 
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Fig. 5. Diagram illustrating "phases" with single-mode and multi-mode operation. 
The data points are cut-off values obtained from analysis of the form indicated in 
Fig. 4. 

4 Conclusion 

We have considered the effective area ^eff of photonic crystal fibers (PCFs) with a 
triangular air-hole lattice in the cladding. Based on extensive numerical studies of PCFs 
with varying air hole size we have constructed a map of the effective area which can be 
scaled to a desired value of the pitch. We have also utilized the concept of effective area to 
calculate the "phase" boundary between the regimes with single-mode and multi-mode 
operation. In this work we have studied PCFs with a triangular air-hole lattice cladding, 
but we emphasize that the approach applies to microstructured fibers in general. 
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